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Abstract This study focuses on liquid water transport in 
wood above the fiber saturation point in the nonhygro- 
scopic region. The liquid water transport of hygroscopic 
porous materials including wood is usually described by 
Darcy’s law. It requires knowledge of capillarity and intrin¬ 
sic and relative permeabilities. In this study, the capillary 
pressure-water relation and relative permeability were 
investigated using experimental data for wood available in 
the literature. The performance of three models (Brooks- 
Corey model, van Genuchten model, and Durner’s bimodal 
pore-size distribution model) was investigated for the capil¬ 
lary pressure-water relation.These models have advantages 
in that each shape parameter has qualitative physical 
meaning for the pore-size distribution. Most species had 
unimodal pore distributions except for aspen, which had 
a bimodal pore distribution. The van Genuchten model 
represented the capillary pressure-water relation better 
than the Brooks-Corey model. Durner’s bimodal model was 
found to be the most appropriate for the capillary pres¬ 
sure-moisture relation of aspen. The relative permeability 
was calculated by using Mualem’s model, which was com¬ 
pared with the value from the Couture model. From the 
results, the liquid water diffusivity divided by intrinsic 
permeability of wood was estimated. This approach may be 
promising for adopting the liquid water diffusivity for the 
numerical simulation of drying and sorption, although there 
might be considerable variation within and between trees. 
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Introduction 

Several models for simulating the transport of heat and 
moisture in wood have been developed recently. Most 
models have been validated by experiment. Among them, 
sophisticated and widely accepted models that have three 
systems of governing equations are for the behavior of 
wood drying processes. 1-3 

Above the fiber saturation point (FSP), wood is non- 
hygroscopic, and most models assume that liquid water 
transport follows Darcy’s law. Bound water movement in 
the hygroscopic range follows Fick’s law. For the simulation 
of liquid water transport in wood, therefore, three material 
indices such as intrinsic and relative permeability as well as 
the capillary pressure-moisture relation (CMR) or the cap¬ 
illary pressure-saturation relation (CSR) are required. For 
the combined bound water movement, three indices such 
as the sorption isotherm, bound water diffusivity, and water 
vapor diffusivity are necessary. Therefore, the simulation of 
moisture transport in the whole moisture range requires six 
material indices that should be measured in independent 
experiments. 

The successful application of heat and moisture trans¬ 
port models depends on the reliability of information on the 
material parameters necessary for the governing equations. 
However, most models contain at least one adjustable 
parameter, which must be determined by independent 
experiment, and so are not completely predictive. In some 
cases, some properties are adopted arbitrarily from other 
studies without investigating the whether the inclusion of 
such data is appropriate. To review and develop heat and 
moisture transport models, it should be noted that sophis¬ 
ticated models do not necessarily provide better prediction 
without minimizing uncertainties in input parameters. 

Heat and moisture transport coefficients are dependent 
on grain orientation, hysteresis, and the local values of 
moisture content and temperature, as well as physicochemi¬ 
cal properties of wood species. Furthermore, the material 
properties of wood are highly variable even within one tree. 
Thus, the information on heat and moisture transport 
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coefficients is very limited because considerable experimen¬ 
tation is required to ascertain them for all possible condi¬ 
tions. Indeed, the measurement of coefficients at a given 
temperature and moisture content is a very time-consuming 
process. 4 This is a major limitation for practical application 
of the models. For liquid water transport in wood, above 
all, it is difficult to find studies that deal with all three prop¬ 
erties mentioned above. 

The purpose of this study was to investigate an appropri¬ 
ate model for the CSR and the relative permeability of 
wood using experimental data available in the literature 
and to estimate the liquid water diffusivity as a function of 
saturation. 


Theory 

The hygroscopic porous material consists of the solid matrix, 
dry air, and moisture. Moisture is the sum of bound water, 
water vapor, and liquid water. In wood science, liquid water 
is called free water. Moisture transport models for isother¬ 
mal conditions can be divided into three types according 
to the driving potential (Table 1). Using chain rules, any 
assumed driving potential can be changed to the equivalent 
potential of moisture gradient. 

In this study, it is assumed that 

1. wood density does not change with moisture content 
(neglecting dimensional change) 

2. gas pressure in wood is constant (neglecting bulk flow of 
fluid) 

3. temperature is constant (neglecting moisture transport 
due to temperature gradient) 

Model 1 


liquid water density, k IW is the relative water permeability, 
K w is intrinsic water permeability, v w is dynamic viscosity of 
liquid water, P c is capillary pressure, D, is diffusivity of 
water vapor, D B is diffusivity of bound water. 

Bound water and water vapor transport in wood are 
interdependent and can be combined, which results in com¬ 
bined bound water and water vapor diffusivity. 5 ' 6 Equation 
1 can be expressed in terms of a single driving potential of 
moisture gradient by using chain rules if the dependence of 
capillary pressure on saturation or moisture content can be 
identified. 

^- = V{(D w + D^D~ B )Vm] (2) 

where D w is liquid water diffusivity and can be expressed 
as 

Pw k IW K w <)I\ dS . 

w p 0 p w dS dm ( ’ 

The relationship between saturation and moisture 
content should be given for converting the CSR into CMR. 
When the moisture content is above the FSP, the water 
saturation (.S') can be calculated by 

r> ^77 /77pcp y, / a\ 

S = -Vm > m FS P (4) 

m max ~ m FSF 

where m FSP is the FSP of wood, which is usually assumed to 
be 0.3. 

For hygroscopic materials such as wood, it should be 
noted that by adopting an irreducible saturation concept 
(m ln > m FS P ), the set of equations may produce a mathemati¬ 
cal singularity as Couture et al. 3 indicated. 

The water dynamic viscosity is calculated by 

f 1540 S 

v, = p.e*p(-m43+j-——-J (5) 


Darcy’s law has been widely used for liquid water transport 
while Fick’s law has been used for bound water and 
water vapor transport. Assuming the driving potentials of 
capillary pressure and moisture gradient, the multiphase 
moisture transport equation for wood can be expressed as 
Eq. 1. 


Therefore, three material parameters ( dPJdm , k IW , K w ) 
should be determined for the prediction of liquid water 
diffusivity under isothermal conditions. 

Model 2 


^ = V . 

dt 


Pw ^ r\v ^ w 


VP c + (D, + D B )Vm 


( 1 ) 


Po V w 

where m is moisture content based on the moisture fraction 
per unit dry wood, p 0 is the oven-dry wood density, p w is the 


In the case of using water potential as the driving 
potential, 




( 6 ) 


Table 1 . Material properties to be measured and estimated by models 


Model 

Driving potential 

Material properties 



References 



Nonhygroscopic region 

Hygroscopic region 

Estimation by inverse method 


1 

Capillary pressure, 
moisture gradient 

CSR. intrinsic and 
relative permeability 

Sorption isotherm, bound water 
and water vapor diffusivity 

- 

1,2 

2 

Water potential 

CSR, sorption isotherm 


Effective moisture 
conductivity 

9,10,18 

3 

Moisture gradient 

- 


Effective moisture diffusivity 

25-27 


CSR, Capillary pressure-saturation relation 
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where y/ is the water potential and K is the effective mois¬ 
ture conductivity. 

The water potential is the energy that is needed to trans¬ 
port unit moisture, which is the sum of the matric and 
osmotic potential. It is usually applied in soil science. The 
relationship between the capillary pressure and the water 
potential is 


_ 2a cos 6 

Pc = -= ~Pw¥ 


( 7 ) 


where cr is the surface tension of water, r is the radius of 
curvature of the air-water interface into the capillaries, and 
0 is the contact angle between the water and the surface of 
the capillary. 

Similarly to Eq. 3, Eq. 6 can be expressed as Eq. 8 from 
the relation between moisture content and water potential 
or capillary pressure to be determined 



-^-Vm 
p 0 dm 



K dP c 
PoPw dm 



( 8 ) 


Comparing Eq. 8 with Eq. 3, the effective moisture conduc¬ 
tivity has the physical meaning involving terms of intrinsic 
and relative permeability. 

k K 

jy _2 ' v rw XY w /n\ 

K = p w - (9) 

Pw 

K cannot be determined directly and so Cloutier and Fortin 7 
determined it using the instantaneous moisture profile 
method in non-steady-state conditions during drying. 


Model 3 

In building physics, a nonlinear diffusion equation is often 
used for simplicity: 

^ = V-(DVm) ( 10 ) 

dt 

where D is the effective moisture diffusivity. 

D involves the combined effects of liquid water, water 
vapor, and bound water diffusivities, and is often expressed 
as an exponential function. It is usually obtained by the 
Boltzmann transformation method during sorption and by 


the moisture gradient method during drying. s The Boltzmann 
transformation method cannot be used for drying because 
it requires the moisture content at one end to be maintained 
at the initial moisture content. 

Recently, Kang et al. 5,6 studied bound water and water 
vapor transport in the hygroscopic region. The present 
study focuses mainly on liquid water transport above the 
FSP. From Eqs. 2, 8, and 10, the effective liquid water dif¬ 
fusivity should be the same in terms of moisture content 
regardless of which driving potential is used. 

D Pw ^rw^w dP c dS _ K dP c ^ 

Po Pw dS dm p 0 p w dm 

It should be noted that dPJdm is negative because capillary 
pressure increases with decreasing moisture content. The 
CSR in the hygroscopic range (below FSP) may be calcu¬ 
lated by the Kelvin equation using the sorption isotherm. 9 


Results and discussion 

The CSR and relative permeability, apart from intrinsic 
permeability, were investigated among three material 
indices for the simulation of liquid water transport. 


Capillary pressure-saturation relation 

The CSR that is usually referred to as the water retention 
curve in soil science cannot be determined theoretically. 
Measurements of the CSR of wood are very scarce when 
compared with studies of other building materials such as 
concrete and soil. Most CSRs of wood that have been mea¬ 
sured are for desorption. Very rarely is the CSR curve 
measured during absorption and the only example may be 
that of western hemlock by Fortin. 10 

The CSR is determined by direct methods (centrifugal 
method, pressure plate and pressure membrane method) or 
indirect methods (mercury intrusion porosimetry, image 
analysis). 8,11 

Table 2 shows empirical equations for the CSR of wood 
during desorption or drying. Spolek and Plumb 4 used models 
similar to those suggested by Brooks and Corey 12 for south- 


Table 2. Empirical equations for capillary pressure-saturation relation of wood 

Equation Constants Species Reference 


P c {S) = a 1 S~ ai 

P C (S, T ) = 101325[fljSexp(-a 2 5) 

+ a, (1 - S)S~“ a ](1 - 2.79-10" 3 (r - 273)) 


PAS,T,p) = a(T) 


«3 + «4Po 

a 5 -S 


+ a 6pQ (1 V + a l + a SpO . 

1 + flioJ 


P c (S,T) = a(T) 



J(S) = 


‘h 

S-a 2 


+ 


°3 

S + a 4 


- a 5 S + a 6 


Oj = 8.4 x 10 4 , a 2 = 0.63 Southern pine 4 

Heartwood a t = 12.12, a 2 = 5.939, a 3 = 0.046, a 4 = 3.7 Pine 13 

Sapwood flj = 1.937, a 2 = 3.785, a 3 = 0.093, n 4 = 1.4 

a 2 = 3150, a 2 = 1.0 x 10 -4 , a, = 1047, « 4 = 3.368, a 5 = 1.02, Softwood 14 

a 6 = 149.8, a, = 52 350, a 8 = 168.4, a 9 = 3150, 

« 10 = 1.0 x 10" 

rij = 9.547 x 10“ 3 , a 2 = 1.028, a 3 = 2.0023 x 10“ 2 , a 4 = 0.0, Bed of glass 15 

a 5 = 1.2 x 10“ 5 , a 6 = 4.415 x 10“ 3 spheres 


o{T) = (77.5 - 0.185(r- 273.15)) x 10“ 3 
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ern pine in the tangential direction using the centrifuge 
method. Bonneau and Puiggali 13 deduced the capillary pres¬ 
sure curves from data obtained by mercury intrusion poro- 
simetry. Using image analysis, Perre and Turner 14 estimated 
them as a function of local density for softwood. In addition, 
Perre and Moyne 15 adopted the Leverett function to fit 
experimental data for a bed of glass spheres that is fre¬ 
quently used in soil science. 

Capillary pressure is strongly nonlinear with moisture 
content although it was plotted as the logarithm function, 
as shown in Fig. la. Power, exponential, or rational func¬ 
tions were used to fit the experimental data for the CSR 
(Table 2). Each expression for nonlinear curve fitting is 
different among the models, and is not generalized. 
However, a power function type of the CSR based on 
Spolek and Plumb 4 is often adopted for softwood drying 
without much correction of parameter values. 16 ' 17 

Unfortunately, the experimental data are not available 
directly and the CSR have been digitized from graphs that 


(a) 



(c) 



can be found in the literature. Figure lb shows the CSR 
measured by the pressure plate method 9 ’ 10 ’ 18 and Fig. lc,d 
shows the CSR measured by the centrifugal method. 19 As 
expected, capillary pressure changes by the pressure plate 
method can be obtained over a wider range of saturation. 
However, most of the experimental data measured by the 
centrifugal method are in a limited range of saturation from 
0.2 to 0.6. 

Aspen distinctly represents the characteristic bimodal 
CSR curves but other species represent typical unimodal 
curves. The measured capillary pressure ranged from 10 3 to 
10 7 Pa depending on saturation. Western hemlock repre¬ 
sents high hysteresis and the capillary pressure was two 
orders of magnitude less than desorption or drying, as 
shown in Fig. lb. At the same saturation, the capillary pres¬ 
sure of heartwood is expected to be higher than that of 
sapwood due to higher pit aspiration rate. On the contrary, 
however, the capillary pressures for the heartwood of 
American elm and cottonwood were less than for sapwood, 


(b) 



(d) 



Fig. la-d. Capillary pressure versus saturation of wood at room temperature and fitted curves of van Genuchten and Durner models, a Empiri¬ 
cal equations, b pressure plate method, c and d centrifugal method. Solid line , heartwood; dotted line , sapwood 
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the cause of which is unclear. It is interesting that the varia¬ 
tion of the CSR is not high compared with intrinsic perme¬ 
ability. The difference in intrinsic permeability is several 
orders of magnitude even within the same species. 20 


Inverse estimation of parameters by CSR models 

The CSR curves can be fitted to an exponential function 
with polynomials. Table 3 shows the best fitted parameters 
for the experimental data (Fig. lb-d), in which saturation 
was converted to moisture content using Eq. 4. Adopting 
the exponential function, however, there are some limita¬ 
tions because no information can be obtained for the pore- 
size distribution and at least three parameters are needed 
for nonlinear fitting. 

Table 4 shows unimodal and multimodal models for the 
CSR. Brooks and Corey (BC) 12 and van Genuchten (VG) 21 
models represent unimodal pore-size distribution models 
having two undetermined parameters. The Durner (DN) 22 
model is often used as a multimodal pore-size distribution 


model, in which the CSR subcurves representing the pore- 
size subsystems can be linearly superimposed by the VG 
model. 

For the BC model, P cb is the air entry capillary pressure 
value and A characterizes the pore-size distribution. Its 
value approaches infinity for a medium with uniform pore 
size, whereas it approaches zero for a medium with a wide 
range of pore sizes. For the VG model, mn or n-q corre¬ 
sponds to A if the product mn remains constant, and a is 
related to the inverse of the air entry value, but strict defini¬ 
tion is not clear. 

Table 5 shows shape parameters for models with the best 
nonlinear curves fitted for the CSRs. The VG models gave 
better fitting results than the BC models in most cases. For 
aspen with a bimodal pore structure, shape parameters of 
the DN bimodal model are given in Table 6. When the CSR 
curve can be fitted with unimodal models, there may be no 
need to use DN multimodal models. However, the Durner 
bimodal model may give better fitting results in less fitted 
cases such as Douglas fir. Using the best fitted shape param¬ 
eters, the fitted CSR curves are shown in Fig. lb-d. 


Table 3. Exponential equation for capillary pressure-saturation relation, P c = expfflj + a 2 m + a 3 m 2 + a 4 m 3 + a 5 m 4 ) 


Method 

Species 

«i 

@2 

°3 

«4 

«5 

R 2 

Pressure plate method 

Aspen 

24.492 

-50.443 

86.929 

-64.735 

16.224 

0.979 


Red pine 

Western hemlock 

16.348 

-0.3909 

-17.761 

21.228 

-7.0789 

0.990 


- 

25.608 

-25.901 

18.418 

-4.748 

- 

0.982 

Centrifugal method 

Absorption 

American elm 

20.097 

-34.215 

32.214 

-8.9183 


0.983 


Sapwood 

17.065 

-2.219 

-3.229 

- 

- 

0.955 


Heartwood 

Baldcypress 

35.385 

-46.594 

22.267 

— 

— 

0.943 


Sapwood 

17.296 

-11.672 

5.338 

- 

- 

0.954 


Heartwood 

Cottonwood 

19.559 

-20.837 

12.596 

— 

— 

0.979 


Sapwood 

15.463 

-1.427 

-0.539 

- 

- 

0.995 


Heartwood 

Sweetgum 

16.862 

-6.161 

1.232 

— 

— 

0.982 


Sapwood 

Yellow poplar 

22.749 

-19.869 

7.658 

— 

— 

0.948 


Heartwood 

Douglas fir 

29.344 

12.998 

-30.596 

6.223 

9.956 

-7.648 



0.984 

0.928 


Table 4. Unimodal and multimodal models for capillary pressure-saturation relation, and its relative permeability by the Mualem model 


Model 

Capillary pressure 

Relative permeability 

Brooks and Corey 12 


k = S p+r+q/x 


U„J 




(Pc > Pcb) 



Mualem 28 

S = [l + (aP c yJ m 


)-(i-s 1 Tj 


(a > 0, n > 1, m = 1 - q/n) 



van Genuchten 21 

k 


( k y 


S = Y,w i S, 


Y^Mi-s^r] 


1=1 

5, = [l + (a i P c rJ"' 

k m = S” 

i=i 

k 

Y, w ‘ a ‘ 


a, > 0, n, > 1, m, = 1 - qln 


v ,=i / 


0 < Wj < 1, Yw,. = l 
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Table 5. Shape parameters for capillary pressure-water relation by the Brooks and Corey and van Genuchten models 


Species 

Brooks and Corey 



van Genuchten 




^cb 

A 

R 2 

a 

n 

R 2 

Softwood 

Red pine 

2000 

0.1531 

0.6978 

7.169 x 10“ 6 

1.921 

0.9848 

Western hemlock 

- 

23000 

0.2648 

0.9176 

9.145 x 10“ 6 

1.464 

0.9896 

Absorption 

164 

0.2397 

0.7765 

2.669 x 10 -4 

2.275 

0.9997 

Baldcypress 

Sapwood 

20490 

0.7873 

0.9842 

4.212 x lO -3 

1.826 

0.9844 

Heartwood 

31496 

0.8832 

0.9934 

2.510 x 10“ 5 

1.969 

0.9945 

Douglas fir 

626 

0.1612 

0.7461 

1.200 x 10“ 3 

1.168 

0.7518 

Hardwood 

Aspen 

1950 

0.2154 

0.8676 

1.740 x 10 -4 

1.282 

0.8854 

American elm 

Sapwood 

2421 

0.1993 

0.9293 

3.906 x 10' 4 

1.201 

0.9300 

Heartwood 

1773 

0.2157 

0.9493 

5.544 x 10' 4 

1.216 

0.9488 

Cottonwood 

Sapwood 

4250 

0.1979 

0.8849 

3.186 x 10“ 5 

1.336 

0.9615 

Heartwood 

5792 

0.2998 

0.9790 

1.604 x 10' 4 

1.304 

0.9796 

Sweetgum 

Sapwood 

2517 

0.1983 

0.9461 

3.756 x 10' 4 

1.200 

0.9414 

Yellow poplar 

Heartwood 

28 

0.1177 

0.9882 

3.380 x 10“ 2 

1.118 

0.9879 


Table 6. Shape parameters for capillary pressure-water relation by Durner bimodal model 

Species 


a x n x 

a? 

*2 

R 2 

Aspen 

0.5966 

2.192 x 10“ 6 2.669 

2.274 x 10 -4 

2.687 

0.9994 

Douglas fir 

0.6273 

2.998 x 10“ 6 1.755 

3.739 x 10" 4 

5.208 

0.9927 


The shape parameter X may change from 0.19 to 11.67 
depending on soils 20 but it ranges from 0.1177 to 0.8832 in 
cases of the wood investigated in this study. According 
to the BC model, yellow poplar might have a wider pore- 
size distribution than baldcypress. Comparing the X of 
the BC model and n of the VG model, X is close to n —1 in 
most wood species. In the case of absorption of western 
hemlock, however, there is a large difference in q of desorp¬ 
tion because q is close to 2 for absorption but 1 for 
desorption. 


Relative permeability 

For the relative permeability of wood, an excellent review 
is provided by Couture et al. 3 Due to difficulties in measur¬ 
ing relative permeability, many methods for estimating it as 
a function of saturation or capillary pressure have been 
developed in soil science and petroleum engineering. It can 
be divided into two groups: those based on the CSR 122122 
and those expressed empirically in a power function of satu¬ 
ration. According to Couture et al., 3 the accuracy of the 
relative permeability determined experimentally remains 
unclear for most materials and arbitrary power functions 
are often used in heat and moisture transport modeling. 
Especially in the case of wood, the experimental data are 
extremely limited except for Tesoro et al. 23 The estimation 


of relative permeability is still open to discussion because 
there is little experimental evidence. 

First, the arbitrary power function is frequently used for 
wood drying 

fc rw =5 B (n > 1,0 < 5 < 1) (12) 

where n is different in the structural direction due to the 
anisotropy of wood (n = 3 in transverse direction and n = 8 
in the fiber direction). 2 

Second, the relative permeability can be estimated using 
the general expression of the Mualem model if the CSR is 
known in advance. 24 


K,=s p 


r rs y 


(13) 


in which the parameters p (<0), q (>0), and r (>1) determine 
the shape of the relative permeability function./? is a lumped 
parameter that accounts for pore connectivity and tortuos¬ 
ity. Therefore, using the generalized Mualem model (Eq. 
13), the relative permeability of the BC, VG, and DN 
models can be obtained as shown in Table 4. In this study, 
we adopt the Mualem model (p = 0.5, q = 1.0, r = 2). 

Third, the partial derivative of capillary pressure with 
respect to saturation is decreasing without bound (~°°) 
because capillary pressure increases indefinitely as satura- 
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tion tends toward zero. This results in indeterminate value 
of liquid velocity. Therefore, the relative permeability of 
the liquid phase was given by using the following form: 3 




min 

S<s< 1 


dS 



with 0 < p < 1 


(14) 


Comparing the predicted drying results with the experi¬ 
mental data for pine, Couture et al. 3 used the same capillary 
pressure-water relation given by Bonneau and Puiggali 13 
and obtained the best agreement by using p = 0.3. 

For the anisotropic materials such as wood, the relative 
permeability may not be dependent on the structural 
direction because the CSR is not dependent on the struc¬ 
tural direction of test samples. 19 Depending on the struc¬ 
tural direction of wood, therefore, the main difference in 
liquid water diffusivity might be due to the intrinsic 
permeability. 

Using the VG model for the CSR, the relative permea¬ 
bility by Mualem and Couture models can be calculated as 
shown in Fig. 2. The nonlinearity of the relative permeabil¬ 
ity with saturation increases as n decreases and r increases. 
Compared with the Mulaem model, the Couture model 
represents higher nonlinearity with saturation. 


Liquid water diffusivity in wood 

To estimate the liquid diffusivity, material properties such 
as the intrinsic and relative permeability, the CSR, and 
density shoud be given (Eq. 11). Flowever, the intrinsic 
permeability of wood varies widely among different speci¬ 
mens of the same species as well as between different 
species. According to Bao et al., 20 the permeability of 
wood is not affected by its changing density between 
and within species. The longitudinal air permeability at 
moisture content 0.1 ranged between 7.49 x 10" 13 m 2 and 
1.85 x 10“ 12 m 2 for softwood, and between 1.80 x 10" 15 m 2 and 
1.33 x 10" 11 m 2 for hardwood. 



Fig. 2. Comparison of estimated relative permeability of wood by the 
Mualem and Couture models (p = 0.5, q = 1). S 3 and S 8 curves refer to 
Eq. 12 


To compare the liquid diffusivity of wood, Eq. 11 can be 
rearranged as 

I\ L:= _l\ L k nL dP L dS_ 
p 0 p w dS dm 

For the VG model, the gradient of capillary pressure with 
respect to saturation is expressed as 


dS 


^£-1 /m _ £-(l +m)/m 


mna 


(16) 


Similary for the Durner model. 


dS 

dP c 


X ~ w i m i n i( a i p c )"[! + («i^c)"] <m ‘ +1> / P c 


(17) 


The maximum possible moisture content of wood is pre¬ 
dicted theoretically if the cell cavities remain constant in 
size during moisture sorption and m FSP is known. 


m„„ = m, 



(18) 


where the fiber saturation point of wood m FSP is usually 
assumed to be 0.3, 0 is porosity, and G„ is dry specific 
gravity. 

The porosity of wood may be calculated from the wood 
and cell wall density using Eq. 4, assuming that the density 
of bound water is equal to the water density, the oven-dry 
cell wall density is 1540 kg/m 3 , and the pore volume remains 
constant with the moisture changes. 


0 = 1 -J^ =1 - 1QOQ x G ° (19) 

r 1540 1540 

Flowever, it is difficult to measure the dry volume because 
the moisture content can change during measuring. There¬ 
fore, the dry specific gravity can be obtained using the moist 
specific gravity. 


G n 


1 -mG. 


( 20 ) 


where G m is the specific gravity of wood based on dry weight 
and moist volume at the given moisture content. The air-dry 
specific gravity at moisture content 0.12 and green moisture 
content of heartwood and sapwood are given in Table 7. 

Figure 3 shows the estimated D W IK„ for Eq. 15. There is 
high variation between wood species. Aspen shows an 
inflection point around saturation 0.6 due to bimodal pore 
distributions. DJK W of softwood is higher than that of hard¬ 
wood in most cases. The liquid diffusivities of red pine and 
cottonwood sapwood were the highest and Douglas fir and 
yellow poplar showed the lowest values for softwood and 
hardwood, respectively. For western hemlock, the average 
liquid diffusivity for desorption was about 7.7 times higher 
than that for absorption. 

There was no significant difference in DJK„ between 
sapwood and heartwood in the case of American elm. For 
cottonwood, DJK W for sapwood was higher than for heart- 
wood. For baldcypress, however, DJK W for sapwood was 
lower than for heartwood. 
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Table 7. Physical properties of wood species 


Species 

Wood handbook 


Dry 

density (kg/nr 3 ) 

Porosity 

Maximum MC 

dS/dm 

Air-dry specific gravity 

Green MC (heart/sap) 

Aspen 

0.39 

0.95/1.13 

409 

0.73 

2.09 

0.56 

Red pine 

0.46 

0.32/1.34 

487 

0.68 

1.70 

0.71 

Western hemlock 

0.45 

0.85/1.70 

476 

0.69 

1.75 

0.69 

American elm 

0.50 

0.95/0.92 

532 

0.65 

1.53 

0.81 

Baldcypress 

0.46 

1.21/1.71 

487 

0.68 

1.70 

0.71 

Cottonwood 

0.35 

1.62//1.46 

365 

0.76 

2.39 

0.48 

Sweetgum 

0.52 

0.79/1.37 

555 

0.64 

1.45 

0.87 

Yellow poplar 

0.42 

0.83/1.06 

442 

0.71 

1.91 

0.62 

Douglas fir 

0.48 

0.37/1.15 

509 

0.67 

1.61 

0.76 


Data from Forest Products Laboratory 29 
MC, Moisture content 


(a) softwood (b) hardwood 




Fig. 3a,b. Estimated liquid diffusivity of wood by the van Genuchten and Durner models (p = 0.5, q = 1, r = 2). D, Desorption; A, absorption; 
S, sapwood; H , heartwood 


Conclusions 

Assuming that the liquid water transport in wood can be 
described by Darcy's law, knowledge of the capillary pres¬ 
sure-saturation relation or the capillary pressure-moisture 
relation, and intrinsic and relative permeabilities is required 
to simulate the moisture transport model. Investigating the 
experimental data of wood available in the literature, most 
species have unimodal pore distributions except for aspen, 
which has a bimodal pore distribution. The van Genuchten 
model represented the capillary pressure-water relation 
better than the Brooks-Corey model. Durner’s bimodal 
model was found to be the most appropriate for the capil¬ 
lary pressure-saturation relation for aspen. For less fitted 
cases such as Douglas fir with unimodal pore distribution, 
Durner’s bimodal model gave better fitting results than the 
van Genuchten model. The relative permeability was calcu¬ 
lated by using Mualem’s model, and was compared with the 
value obtained from the Couture model. From the results, 
the liquid water diffusivity divided by intrinsic permeability 
of wood can be estimated as a function of saturation or 


moisture content. This approach may be promising for 
adopting the liquid water diffusivity for numerical simu¬ 
lation of drying and sorption and to minimizing experi¬ 
mentation required for measurement of liquid transport 
coefficients. In the future, more studies are required on the 
effects of physicochemical properties of wood on liquid 
transport parameters because there may be considerable 
variation within and between trees. 
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